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ABSTRACT 

We present resolved stellar photometry of 4 fields along the major axis of the M33 disk from images 
taken with the Advanced Camera for Surveys aboard the Hubble Space Telescope. The photometry 
provides a detailed census of the red clump in all fields and reaches the ancient main sequence in 
the outermost field. Through detailed modeling of the color-magnitude diagrams, we find that the 
percentage of the stellar mass formed prior to z = 1 changes from 71±9% in the inner-most field 
to 16±6% in the outermost field. The disk shows a clear trend of increasing scale-length with time, 
evolving from = 1.0 ±0.1 kpc 10 Gyr ago to = 1.8 ±0.1 kpc at times more recent than 5 Gyr ago, 
in agreement with analytical predictions for disk growth. Beyond the disk truncation radius, however, 
the stellar density profile steepens with time and the age gradient reverses, in agreement with recent 
simulations. The late and slow growth of the stellar disk may be due in part to the low mass of M33. 
Subject headings: galaxies: individual (M33) — galaxies: spiral — galaxies: evolution — galaxies: 
stellar content 



1. INTRODUCTION 

The scale-lengths of galaxy disks are thought to in- 
crease with time due to late time accretion of gas at 
large radii and exhaustion of gas in the centers of galaxies 
(lLarson|ll 97"ilMatteucci fc Francoi3ll989tlBurkert et all 
[1992; Chi appini et al.lll997t iNaab fc Ostrikeill2006l and 
many others). This process, which has been called 
"insi de-out" growth repro duces many of the physical 
(e.g.. IBouwens et"aI]|1997D and chemical properties of 
disks. Observations of color gradients in spiral galax- 
ies suggest that the mean ages of the stars d ecrease from 
the inside-out (e.g.. IMacArt hur et al.l l2004[ ). and obser- 
vations of extended ultra-violet disks show that very re- 
cent star f ormation is prevalen t in the outer regions of 
disks (e.g.. IBoissier et aT1l2008D . While these trends are 
consistent with inside-out growth, broadband colors can 
be biased by variations in star formation, dust extinc- 
tion, and metallicity, while ultra-violet emission probes 
only the current disk properties. A more direct approach 
to detecting the products of long-term inside-out disk 
growth is through measurements of the age distribution 
of the resolved stellar populations. 

The inside-out growth of stellar disks may be most 
easily detected in low-mass disks. If massive disks 
formed earlier than low-mass disk s (i.e. "down-sizing" 
Cowie et al.lll996HLillv et al.ll2003HNeistein et al.ll2006l : 



Fontanot et al mm . and many others), the stars in mas- 



sive disks have had more t ime to be mixed by radial redis- 
tribution processes (e.g.. lRoskar et al.ll2008l ). Therefore 
detecting the signature of inside-out growth using stellar 
populations may be simplified in low-mass disks. 
M33 is the nearest example of a low-mass disk galaxy. 
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It is also a bulge-less system (jMcLean fc LidlT996[ ). so 
that there is no confusion between disk and bulge popu- 
lations. If the stellar disk of M33 grew from the inside- 
out, the age of the stellar populations should vary with 
radius, with a higher percentage of ancient stars in the 
inner disk. Indeed, ground-based observations of M33's 
bright stellar populations suggest that its outer disk may 
have formed at late epochs and may even be on going 
(lDavidgell200l iBlock et aT]l200l iRo^e et al.ll2005l) . On 
the other hand. M33 has a truncated disk (at r~8 kpc, 
iFerguson et al.l 120071 ) and the age of populations of the 
outer disk/halo sho w an age increase w ith radius based 
on HST/ACS data (|Barker et al.ll2007bl lal) . in contrast to 
the expectations of inside-out growth. 

In this letter, we use data from the Hubble Space 
Telescope (HST) Advanced Camera for Surveys (ACS) 
to study the age distribution of the stellar populations 
of the M33 disk as a function of galactocentric radius. 
These distributions show a clear trend: regions nearer 
to the center formed a higher fraction of their stars at 
earlier times leading to a steady growth in the disk scale 
length. Our measurements therefore provide strong ob- 
servational support for inside-out disk growth. 

More rigorous checks of this analysis and more detailed 
studies of these fields will be provided in a future full- 
length paper (Holtzman et al. 2009, in prepa ration) . We 
assume a distance of 800 kpc (|Lee et al.ll2002l ) for conver- 
sions of angular measuremen ts to physical distances. We 
assume an incHnation of 56° (|Corbelli fc SchneidedfT997h 
for de-projections. We adopt a WMAP^ cosmology for 
all conversions between time and redshift. 



2. DATA ACQUISITION AND PROCESSING 

We obtained HST/ACS imaging of 4 inner fields along 
the major axis of the M33 disk as part of the HST GO 
program 10190 (see left panel of Figure [T|). The obser- 
vations were performed from 09-Sep-2004 to 20-Feb-2005 
through the F606W (wide-F) and F814W (/ equivalent) 
filters. The total exposure time in the innermost field 

^ http: / /lambda. gsfc.nasa.gov/product/map/dr2/params/lcdm_wmap.cfm 
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was 6260 s and 6482 s in F606W and F814W, respec- 
tively. The total exposure time for each of the other 3 
fields was 21140 s and 26300 s in F606W and F814W, 
respectively. 

We also analyzed HST/ACS imaging of 3 outer fields 
along the minor axis of the M33 disk (GO-947 9 ), whic h 
have been previously studied by [Barker et al.l (|2007bD . 
Our reanalysis of these fields allowed us to extend our 
radial profiles to larger radii and provided a consistency 
check between our analysis technique and those used for 
previous work. 

Resolved stellar photometry was measured using the 
photometric pipeline o f the ACS Ne arby Galaxies Sur- 
vey Treasury (Dalcanto n et al.ll2009r ). The pipeline uses 
the DOLPHOT (Dolphin ' 2000") software package to mea- 
sure point spread function photometry and to deter- 
mine errors and completeness as a func tion of color and 
magnitude using artificial star tests (jPalcanton et al.l 
I2009t iWiUiams et"all l2009f) . The final stellar catalogs 
contained 414780, 386894, 407741, and 309249 stars, 
from the innermost field to the outermost field, respec- 
tively. The resulting color-magnitude diagrams (CMDs) 
are shown in Figu re [21 Our measurements of the 
[Barker et all (|2007bf) fields yielded 13522, 3958, and 1548 
stars from the innermost field to the outermost field. 
The resulting CMDs were equivalent to those shown in 
[Barker et alT(|2007br ). 

The CMDs shown in Figure [D qualitatively show the 
relative number of main sequence stars to giants increas- 
ing with the distance from the center of M33, suggest- 
ing decreasing stellar ages with increasing radius. To 
quantify this trend, we measured the star formation 
history (SFH) of each of the 4 field s using the soft- 
ware package MATCH (|Dolphinll200l l. which finds the 
best-fitting distribution of stellar ages and metallicities 
that reproduces the observed CMD assuming a Salpetcr 
(19551) initial mass f unction (IMF) and t he stellar evo- 
lution isochro nes of [Girardi et all ()2002f ) including re- 
cent updates (jMarigo et all [2008D . The details of the 
fitting me thod used we re the same as those employed 
in Willia ms et al.[ ()2009f ). where results at old ages were 
binned to low time resolution to reduce errors and avoid 
over-interpreting details that could be due to deficiencies 
in the stellar evolution models. This technique for recov- 
ering the SFH from the CMD is well-tested and is found 
to be robust against differences in stellar evolution mod- 
els and fitting techniques when the distance and mean 
extinction to the field are allowed to be fitted as free 
parameters. 

Uncertainties were determined by fitting the CMDs 
with a range of assumed distance and mean reddening 
values. These uncertainties were then added in quadra- 
ture to the standard deviation calculated from Monte 
Carlo tests where the observed CMD was resampled and 
refitted one hundred times. 

A constant mean extinction was applied to all 
model stars, u sing t he ACS extinction coefficients from 
[Girardi et al.l ()2008l ). Greater reddening was assumed 
for stars with ages <100 Myr, but these ages are not 
relevant to this study. As a first-order attempt to ac- 
count for differential reddening, we applied a uniform 
spread in Ay to the model photometry of 0.6, 0.4, 0.2, 
and 0.0 magnitudes for the innermost to outermost fields 
(respectively). These values improved the quality of the 



model fits compared to fits with a single uniform redden- 
ing value for each field. However, differential reddening 
did not significantly alter the inferred SFHs (see § 3.4), 
which are shown in Figure [3] 

3. RESULTS 

3.1. Star Formation Histories 

The M33 disk is naturally divided into two compo- 
nents: the inner disk and outer disk. These components 
are divided by a break in the ex ponential density pro - 
file located at a radius of ~8 kpc (|Ferguson et al.ll2007f) . 
The radius of this truncation point is located between our 
outermost field and the innermost [Barker et al.l ()2007b[ ) 
field. 

The SFHs of the fields (see Figure^ show that the bulk 
of the stars in each ACS field formed at different times, 
as can be seen most easily by looking at the cumulative 
fraction of stars formed during each time bin. This value 
is plotted as a function of age for all 4 inner fields in 
the right panel of Figure [TJ The plot clearly shows that, 
while the majority of the stars near the center of the disk 
had formed by 2; = 1, the bulk of the stars farther out in 
the disk formed later. 

The growth is also shown in Figure HI which plots the 
reconstructed stellar mass density profile of the M33 disk 
at 5 different epochs. This plot shows that most of the 
stellar disk growth occurred at radii of ~3-7 kpc, rather 
than at the center, leading to an increase in the scale 
length of the inner disk (r<8 kpc) with little change in 
the central stellar surface density. In contrast, while the 
density profile of the inner disk has fiattened as the disk 
has aged, that of the outer disk has steepened. 

3.2. Implications for the Inner Disk 

We have fit the inner disk (r<7 kpc) with an exponen- 
tial profile and plotted the results in the central panel of 
FigureHl The scale length of M33's disk has increased by 
a factor of ~2 from 10 Gyr ago to the recent past. Most 
of this stellar disk growth occurred from 10 Gy r to 5 Gyr 
ago. This growth is similar to the predictions of [Mo et al.l 
(1998), which are based on simple scaling models for the 
evolution of the virial mass and virial radius. 

The current stellar m ass scale- length is lo nger than 
that of the /sT-band hght (jRegan fc Voge]|[1994 5.8' in K 
vs. 7.7' here), suggesting a changing M/Lk with radius. 
This discrepancy has been seen in a dynamical es timate 
of the scale length of M33's mass (|Ciardullo et al.ll'2004j) . 
However, the i^T-band scale length is similar to the scale 
length we calculate for stellar populations older than 5 
Gyr. 

The right panel of Figure |4] shows the extrapolated 
central surface density of the exponential fits to the stel- 
lar mass density profile. The uncertainties are signifi- 
cant, and the data are consistent with no evolution in 
the central surface density of the stellar disk (solid line). 
We compare the central surface d ensities w i th the expec- 
tations from the scaling laws of [Mo et all (|1998f ): how- 
ever, those laws track only the mean distribution of the 
baryons and formally predict a decrease in surface den- 
sity with time. This prediction is not easily compared 
with our measurement for the stellar disk, for which such 
a decrease is likely to be unphysical. 

Based on Figure 31 we conclude that most of the 
stars outside of the innermost field formed since z — 
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1. However, this conclusion differs from what has 
been seen in the outer disk of more niassive spirals, 
such as M31 (jBrown et al.ll2006[ ). M81 (jWilliams et all 
I2009D . and the Milky Way thick disk, in which most 
of the stars appear to have formed prior to z = 1. 
This early formation of the disks of massive spirals 
is also seen in redshift surveys, which ar e sensitive 
to equally massive spirals beyon d z = 1 (iLilly et al.l 
1998; Ravindranath et al. 200 41 iPapovich et alJ 120051 : 
Sargent et al.l I2007t IMelbourne et al.ll2007[ ). Such sur- 
veys are not sensitive to low- mass disks like M33. The 
difference between our results for M33 and more massive 
disks may be due to the fact that such low-mass galax- 
ies are thought to perform most of their star formation 
at later epochs, an effect that has been deemed "down- 
sizing." 

3.3. Comparison to the Outer Disk 

There is a striking difference between the age gradi- 
ent seen in our data for the inner disk and that seen in 
the far outer disk bv I Barker et al.l ([2007b). Within their 
series of outer disk fields, their measurements show the 
mean age of the stars increasing with radius. In con- 
trast, our analysis shows the opposite effect within the 
break radius, marked by the vertical gray line in the left 
panel of Figure SI Similar differences between the inner 
and outer disk are seen in current simulations of disk for- 
mation, where radial mixing redistributes older stars to 
larger radii. This mixing produces a break in the expo- 
nential disk profile, inside of which stellar ages decrease 
with increasing radius and outs ide of which stellar ages 
increase with increasing radius (jRoskar et al.ll2008f l. Al- 
though this radial mixing scatters stars into the outer 
disk, the majority of the stars in the inner disk remain 
close to their formation radius. Thus, while mixing pro- 
duces the break and age gradient in the outer regions, 
it is not sufficient to erase the radial age gradient estab- 
lished during the formation of the inner disk. 

Our independent reduction of the I Barker et al.l ()2007bl ) 
data and derivation of the SFH reproduces their observed 
increase in mean age with radius in the outer disk. The 
resulting surface densities of these outer fields are in- 
cluded in Figure ID Fits of a single exponential to the 
full broken exponential (see center panel of Figure S]) av- 
erage the inner and outer disk gradient and thus show no 
significant change in scale-length with age. This exercise 
demonstrates the importance of measuring profiles of the 
disk inside and outside of the disk break separately. 

Finally, our results show some very old stars at all 
radii in the disk. However, in disk formation and evolu- 
tion simulations, radial redistribution places stars of the 
oldest ages throughout the disk, in cluding much large r 
radii than where they were formed (jRoskar et al.ll2008f ). 
suggesting that the oldest stars seen at larger radii could 
have been formed closer to the galaxy center. The old- 
est stars may also have been accreted. In any case, the 
inner disk apparently had a shorter scale-length in the 
past, and the trends seen in mean age and scale- length 
in our measurements mimic those seen in simulations and 
semi-analytic models of inside-out growth. 

3.4. Reliability Tests 

Because our conclusions rely strongly on the reliability 
of our technique to measure SFHs from resolved stellar 



photometry, we have performed several tests to check 
the robustness of our measurements. Our conclusions 
draw upon the age distribution of stars with ages of 2-12 
Gyr, resolved into 4 time bins. We performed multiple 
binnings to check the sensitivity of our results to bin- 
ning scheme and verified that the same age gradient was 
present even for larger more conservative time bins than 
those finally adopted. 

While such high time resolution in the 2-12 Gyr range 
is more robust at the depths of the outer 3 major axis 
fields, the ability to distinguish between stars with ages 
of 5 and 10 Gyr may be less reliable in the shallower in- 
nermost field. We therefore calculated the scale length of 
the inner disk as a function of time excluding the inner- 
most field and found the trend to be equivalent to that 
measured when the innermost field is included. 

We also tested the effects of forcing the distance to 
be the same in each fit by assigning a fixed distance of 
(m — M)o=24.6 to all of the fields. This test decreased 
the scale lengths at all epochs, but left the trend the 
same. 

The CMDs contain stars with a spread in reddening, 
such that the red clump appears slightly smeared along 
the reddening line. To test the effects of our assumption 
that differential reddening effects are less severe in fields 
farther out in the disk (see § 2), we refit our data first 
assuming the same spread in reddening (Ay = 0.6) in 
the innermost 3 fields and again assuming no differen- 
tial reddening in any of the fields. While these changes 
affected the resulting overall scale lengths, causing the 
trend to be from l.OiO.l kpc to 1.7±0.2 kpc for a spread 
of Av = 0.6 and from l.liO.l kpc to 2.0±0.2 kpc for no 
spread at all, the changes did not significantly alter the 
resulting SFHs or change the trend of increasing scale 
length with time. More sophisticated tests of these ef- 
fects will be reported in the more comprehensive paper 
on this dataset (Holtzman et al. 2009, in preparation). 

We also checked for consi stency with previous work. 
Our measured SFHs for the [Barker et all (|2007b[ ) fields 
were consistent with theirs within the errors except for 
a factor of ^2 offset in the nor malization, due to the 
IMF. Since [Barker et all ()2007bD assumed an IMF that 
flattens at low masses, their models required lower star 
formation rates to reproduce the observed CMD. When 
normalized to produce the same number of stars from 1 
to 100 M0, the ratio of the total number of stars pro- 
duced by their adopted IMF to that produced by our 
adopted IMF is '^0.5. Correcting for this offset puts our 
resulting SFHs in good alignment with theirs. Changes 
in the IMF and mass cutoff serve only to shift all of the 
surface density data in Figure [4] up or down. The scale- 
lengths and fractional changes in central surface density 
are thus unaffected by t hese IMF issues, un less the true 
IMF varies with radius ijMeurer et al.ll2009f ). 

4. CONCLUSIONS 

We have modeled the CMDs of 7 HST/ACS fields in 
the M33 disk at a range of galactocentric radii. Four 
of these fields are newly-observed and lie along the ma- 
jor axis of the inner disk. The resulting stellar age dis- 
tributions demonstrate that, within the inner disk, the 
age of the majority of the stars decreases with increas- 
ing distance from the galactic center. In addition, most 
of the disk stars outside of ~3 kpc formed since z = 1, 
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later than what is seen in more massive spiral galaxies 
both in stellar populations studies and in redshift sur- 
veys. Comparisons of the stellar population properties in 
th ese fields with those o f the outer M33 disk as measured 
bv [Barker et al.l ()2007bl ) show an inversion of the radial 
dependence of stellar age, providing support for some 
current simulation results. Together, these results pro- 
vide strong observational support to inside-out growth 



of low-mass spiral stellar disks, down-sizing, and the sig- 
nificance of disk breaks in formation models of galaxy 
disks. 
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Fig. 1 . — Left: A digitized sky survey image of M33. Overplotted in color are our 4 HST/ACS fields. The inset image shows the locations 
of the 3 [Barker et al. (2007b) fields as red squares. Right: The cumulative fraction of stars formed in the four inner fields as a function of 
time. The colors of the lines are matched to the colors of the field outlines in left. 
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Fig. 2. — Upper Left: The F606W-F814W CMD of the innermost field. Contour levels indicate density of points where the plot would 
be saturated (levels are 1,1.5,2,2.5,3,4,6,8,12,16,20 xlO* mag~^). Photometric errors as a function of F814W magnitude are shown on the 
left side. The field name and reddening vector are shown in the lower-right corner. Upper Right: Same as upper left, for the next field out 
from the center. Lower Left: Sa me as upper r ight, fo r the next field out. Lower Left: Same as lower left, but for the outermost field. For 
CMDs of the archival fields, see iBarker et al.l ll2007blV 
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Fig. 3. — Upper Left: The star formation rate as a function of time of the innermost field. A thicl<, gray, dashed Hne denotes the mean 
rate for the field. Upper Right: Same as upper left, but for the next field out from the center. Lower Left: Same as upper rig ht, but for the 
next fie ld out. Lower Left: Same as lower left, but for the outermost field. For SFHs of the archival outer fields, please see [Barker et all 
Il2007bl) . 
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Fig. 4. — Left: M33 stellar surface densit ies as a function of de-projected radius as computed from our measured SFHs (stars) and from 
our re-analysis of the [Barker et all Il2007bl') fields (diamonds). Different colors show the total stellar densities that would be observed at 
different lookback times, noted in the corners. The absolute mass density values on the ordinate can shift based on the IMF. The best-fit 
exponential scale-length for each epoch using just the fields inside the disk break is listed in the upper-right. Those using all of the fields 
are listed in the lower-left. Data po ints have been offset f rom one another by 0.1 kpc to avoid overlapping error bars. The vertical gray line 
marks the disk break measured by [Ferguson et aTl 1)20071 ) . Center: The disk scale length as a function of age inside the disk break (solid 
line) and from all fields (dotted line). The dashed line shows the scaling relation of lMo et al.l 11199 81 normalized to our 5 Gyr measurement. 
Right: The central surface density as a function of age as measured from fields inside the disk break (solid line) and from all fields (dotted 
line). The dashed line shows the scaling relation of lMo et al.l l l 19981) normalized to our 5 Gyr measurement. 



